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Abstract 

In the limit of small velocities, the dynamics of half-BPS Yang-Mills-Higgs solitons 
_C ■ can be described by the geodesic approximation. Recently, it has been shown that 

quarter-BPS states require the addition of a potential term to this approximation. We 
explain the logic behind this modification for a larger class of models and then analyse in 



detail the dynamics of two five-dimensional dyonic instantons, using both analytical and 



numerical techniques. Nonzero-modes are shown to play a crucial role in the analysis 
of these systems, and we explain how these modes lead to qualitatively new types of 
dynamics. 
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1 Introduction 

Much about the kinematics of supersymmetric solitons is known by virtue of the fact that their 
properties are related to the superalgebra of the underlying field theory. On the other hand 
the dynamical aspects of these objects are much harder to analyse. Extracting information 
about the dynamics of solitons by direct analysis of exact, time dependent solutions to the 
full equations of motion is generically too complicated. Instead, one uses the approximate 
but very powerful moduli space technique first introduced by Manton [1] for the study of the 
dynamics of four- dimensional Yang-Mills-Higgs magnetic monopoles. 

Manton's method is based on the idea that the dynamics of solitons can be studied, to 
lowest order in the velocity, by a quasi-static approach. The set of static solutions, modulo 
small gauge transformations, is parametrised by a number of collective coordinates or moduli, 
and the space they span is called moduli space. At each moment of time the system is assumed 
to be infinitesimally close to a particular static configuration on this space. An effective 
action which tells us how the system evolves from one static configuration to another can be 
constructed as follows. 

One first promotes the parameters in the static solution to be slowly varying functions 
of time. In order to deal only with the physical degrees of freedom, one has to ensure that 
this configuration, with time- dependent parameters, also satisfies the Gauss' constraint at 
the lowest order in the velocity and at each point in time. Very often this requires adding 
compensating gauge transformation terms to the static solution. This time-dependent con- 
figuration does by construction not change the potential energy, and the deviation from the 
original static solution is therefore called a zero-mode. If the initial (static) configuration was 
minimising the potential energy, then the zero-modes are deformations which will make the 



slowly moving system evolve through the set of states with the same, minimal value of the 
potential energy. 

By inserting now this time-dependent solution into the action, expanding it to lowest order 
in the velocity and integrating over the space-like coordinates, one arrives at an effective action 
for the moduli of the static configuration. This action is just an ordinary non-linear world- 
line sigma model with the moduli space as the target space. The metric on it is naturally 
induced from the field theory action 1 , and the dynamical evolution of the system is described 
by geodesies on the moduli space. Although strictly speaking this approach is valid only in the 
limit of very small velocities, an analysis of the dynamics of two Yang-Mills-Higgs monopoles 
by Manton and Samols [3] shows that this method is surprisingly accurate. Radiation effects 
contribute only a tiny correction at speeds up to a considerable fraction of the speed of light. 

An interesting feature that the moduli space method has brought to attention is the fact 
that four-dimensional 1/2-BPS dyonic solitons can be interpreted as zero-mode excitations of 
an underlying purely magnetically charged soliton. In particular the Julia-Zee dyon [4] arises 
naturally in this approach as an excited 't Hooft-Polyakov monopole. This can be seen by 
expanding the Julia-Zee field theory solution to lowest order in the electric charge. One is 
then left with an unmodified 't Hooft-Polyakov monopole and a non-vanishing A component 
of the gauge field. Since the electric part of the energy of the dyon is a purely "kinetic" 
contribution to the total energy, the question arises as to what it is that is moving in these 
solutions. It turns out that the A component of the solution is precisely equal to the global 
gauge rotation zero-mode of the starting 't Hooft-Polyakov monopole, with the velocity equal 
to the electric charge. The direction of this additional zero-mode is compact, which accounts 
for the quantisation of the electric charge. One concludes that the moduli space approximation 
for monopoles automatically incorporates dyons, and dynamical processes can lead to charge 
exchange. 

Recent developments in the study of dualities in gauge theory and string theory have 
uncovered novel four and five-dimensional 1/4-BPS dyonic solutions, starting with the work 
of Bergman [5]. In contrast to the 1/2-BPS dyons discussed above, the dynamics of these 
solitons is naturally described on the moduli space of approximately static configurations. 
This is strongly tied to the way in which those solutions are constructed in the first place. In 
the cases known so far, the 1/4-BPS dyons can be seen as deformations of 1/2-BPS electrically 
uncharged solutions, obtained by turning on an additional Higgs field. The equations one has 
to solve are an unmodified BPS equation for the underlying 1/2-BPS configuration, together 
with an equation for the additional scalar field in the background of the underlying 1/2-BPS 
soliton. The latter has a unique solution, and thus provides a one-to-one map between the 
moduli spaces of the two types of solitons. However, turning on the additional Higgs makes the 
system unstable, because the Higgs field introduces a new attractive force which is unbalanced. 
By further adding an additional electric charge, a compensating force is introduced which 
prevents the configuration from collapsing. The resulting solution is dyonic in the sense 
that it carries both the charge of the original 1/2-BPS soliton as well as the electric charge 
introduced to prevent collapse. 

The instability of the intermediate step, where the Higgs is turned on but the electric 



1 Remarkably, the metric on the most studied example of two SU(2) monopoles was actually never obtained 
this way. Instead, Atiyah and Hitchin [2] were able to fix the metric completely using only symmetry arguments. 
This method has clear calculational advantages, as the zero-modes for this case are very complicated. In general 
situations, however, there are not enough symmetries and one has to resort to the zero-mode method described 
above. 



charge is still zero, can alternatively be seen through the appearance of a potential on moduli 
space. This potential can in principle be derived just like the metric part of the action, but 
most of what is known is actually based on supersymmetry. Since rather general arguments 
show that the potential has to be written as the square of a Killing vector on moduli space, 
the analysis of Alvarez-Gaume and Freedman [6] can be used to argue that the orbits of the 
Killing vector have to preserve the three complex structures in order for it to be compatible 
with the four supersymmetries. This observation, due to Tong [7] and elaborated on in many 
other papers, has led to a determination of explicit potentials for several dyon systems in 
SU(3) Yang-Mills-Higgs. The general structure of the effective action, and its supersymmetric 
extension, has been derived in a series of papers by Bak et al. [8, 9] and Gauntlett et al. [10, 11]. 
The electric charge again arises from a dynamical effect on the moduli space, just as 
in the case of the Julia-Zee dyon. However, the charge deformation is now a zero-mode of 
an approximately static configuration. Even though the starting configuration is unstable, 
the system becomes stable once the electric charge zero-mode is added. In fact, this excited 
configuration corresponds to the exact 1/4-BPS solution. Heuristically, the stabilisation arises 
because the centrifugal force due to the motion of the system in the compact direction is 
balanced against the force produced by the potential. In summary, these new static 1/4-BPS 
dyons are identified with closed geodesies of constant energy on the moduli space of 1/2-BPS 
solitons with a potential 2 . The conserved quantity specifying different geodesies is the angular 
momentum, and its value is identified with the electric charge carried by the dyon. 

Although the origin and the form of the potential are thus rather well understood, it is 
important to realise that inclusion of a potential leads one beyond the strict moduli space 
approximation. The difference becomes most apparent when one explicitly tries to analyse 
the dynamics of such models, and the present paper is an attempt to do so. The addition of the 
potential clearly only makes sense when its magnitude and slope are small compared to other 
scales in the problem. In this regime, the unstable configurations solve the equations of motion 
to lowest order in the Higgs expectation value, and small deformations of this approximate 
solution can be studied. The picture one should have in mind is very similar to the picture of 
the 1/2-BPS case: the system is restricted to move at the bottom of a potential valley, which 
is surrounded by steep mountains. The key difference is that in the 1/4-BPS case the valley 
is slightly tilted. The smallness of the slope is required to prevent the system from gaining 
enough kinetic energy to start climbing the mountains. 

With such a tilted valley, the deformations one has to consider are not just the zero- 
modes (i.e. the modes which keep the potential energy constant) but also the destabilising 
modes (i.e. the modes which change the potential by moving slowly down the tilted valley). In 
general it is very difficult to include nonzero-modes in the dynamics, because there is no simple 
way to determine which of the nonzero-modes are most relevant. While a systematic approach 
has been described by Manton [12], there are not very many models for which this has been 
worked out in detail (see Manton and Merabet [13] for one application and further references). 
In the present case, however, one can argue that an appropriate set of nonzero-modes consists 
of the lifted zero-modes of the underlying 1/2-BPS soliton. 

Given this insight, one can study the explicit solutions to the equations of motion of 
the effective action, and the symmetries that govern them. We will show how the inclusion of 
nonzero-modes sometimes leads to the appearance of several separated scales in the low-energy 
dynamics, a new feature which does not occur for the 1/2-BPS models. Instead of ignoring 



2 For lack of a better alternative we keep referring to the solutions to the equations of motion as "geodesies" . 
Of course these solutions are not necessarily paths of extremal length when the effect of the potential is included. 



the nonzero-modes, one can actually integrate them out explicitly, and the resulting dynamics 
is qualitatively different from the dynamics predicted by just the zero-modes. Throughout 
the paper we will use the explicit example of two dyonic instantons, which are solitons in 
4+1 dimensional Yang-Mills-Higgs theory. They have the advantage over other 1/4-BPS 
solitons that the metric on moduli space is not too complicated, and can be constructed 
explicitly from a paper by Osborn [14]. 

The plan of the paper is then as follows. First, we recall the details of 1/4-BPS (multi) 
dyonic solitons in 4+1 dimensional Yang-Mills-Higgs theory. In section 2.2 we then exhibit 
the appearance of a potential on moduli space, illustrated by this particular example. The 
dynamics of a single dyonic instanton can be studied analytically, which we do in section 2.3. 
The main part of our paper, contained in section 3, is concerned with the study of the two- 
soliton sector. We first derive the metric and potential, and then perform a detailed analysis 
of the dynamics. Qualitatively new features are found in section 3.3 and the final section is 
concerned with the separation of scales in the low-energy theory. We end with some comments 
on work in progress concerning 1/4-BPS dyons in 3+1 dimensions. 

Two appendices have been added. One discusses technical details of the two-instanton 
moduli space, while the other one summarises our conventions on quaternions and lists a 
number of useful expressions involving them. 

2 Dyonic instantons and moduli spaces 

2.1 The static dyonic instanton solution 

After the general introduction we will, in the rest of the paper, specialise to a concrete example, 
namely the study of dyonic instantons. Before we discuss their dynamics, let us first introduce 
the static solution. The model under consideration is the five dimensional Yang-Mills-Higgs 
action, given by 



S = - fd 5 x (i tr {F^F^) + \ tr (£>„$£>"$) ) , 



T, 



with \x = (0, (i = 1..4)). This model can be extended to a supersymmetric one, but we will 
not need the details of this extension in the present paper. 
For static solutions the equations of motion reduce to 

DjF* + [A , D l A ] - [$, D% = , (2) 

DjF i0 - [$, D°$] = , (3) 

DiD'Q = . (4) 

The first line clearly shows that by choosing A = ±<l> the last two terms cancel and the 
remainder of the equation is solved by an (anti) self-dual field-strength 3 . One can alternatively 
see this by rewriting the Hamiltonian in Bogomol'nyi form, which leads to the following 
BPS equations, 

F 0i = s'D t $ , (6) 



5 For Julia-Zee dyons all three terms mix and the relation between A and $ is less rigid. 



with s and s' two arbitrary signs (corresponding to the sign of the instanton and electric 
charge respectively). Combining the second BPS equation with the Gauss law and going 
to the Aq = ±$ gauge one arrives again at the equation of motion (4) for the Higgs field, 
which is just a Laplace equation in the background of the instanton. Hence as advertised, the 
equations for the 1/4-BPS dyon reduce to two equations: a BPS equation for the instanton 
and a Laplace equation for the scalar. An important observation is that for each instanton 
solution to (5) there is a unique solution to (4). Hence there is a one to one map between the 
moduli space of instantons and the moduli space of dyonic instantons. 

For instantons with identical group embedding angles, the solution to these equations is 
given by Lambert and Tong [15], extended by Eyras et al. [16] and given a brane interpretation 
by one of the authors in [17]. For s = — 1 one obtains, 



(7) 



where the function H is given by 
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The symbols rjfj are the 't Hooft symbols and the matrices T a are Pauli matrices. This solution 
is characterised by two conserved charges: the instanton number 

I=ljd*xe ijkl tr(F ij F kl ) = N, (9) 

and the electric charge 

q = s'- fd 4 x tr( D&) 2 = 4n 2 vs' J^ Pi ■ ( 10 ) 

i 

This solution can be extended to an arbitrary ADHM configuration for the four- dimensional 
subspace, as the solution to the Higgs field in the instanton background has been derived in 
terms of ADHM data by Dorey et al. [18]. We will make use of this later when we study the 
two-dyon dynamics. 

2.2 Modes, zero-modes and the moduli space approximation 

We now want to analyse the dynamics of the dyonic instanton solutions. The solution which 
was described in the previous section is an exact, static solution of the equations of motion. In 
principle one could try to construct the zero-modes of this configuration and the correspond- 
ing effective action for them. This method was followed by Bak and Lee [19] for 1/4-BPS 
monopoles. Our approach is more in spirit of Manton's original treatment [20] of Julia-Zee 
dyons 4 . Namely, we will try to make a link between dyonic instantons and instantonic solitons 



4 It is an interesting open problem to construct a moduli space effective action for two Julia-Zee dyons 
starting from the dyonic solution, and understand to what extent such an action agrees with the one constructed 
with monopoles as starting point. This would be the analogue of starting from two dyonic instantons in the 
model under consideration here. 



(i.e. instanton solutions trivially embedded in the 4+1 dimensional Yang-Mills theory) and 
profit from the vast knowledge about zero- modes of the latter 5 . As we have already mentioned, 
there is a one-to-one map between the moduli space of instantons and the moduli space of 
dyonic instantons. However, since the non-zero Higgs expectation value of a dyonic instanton 
can not be obtained as a (zero) mode deformation of the instanton, one cannot start from the 
instanton moduli space. Instead one takes the "instanton solution" in the Yang-Mills-Higgsed 
theory as a starting point. 

It is well known that due to Derrick's theorem [21], turning on the scalar field in the 
four-dimensional Euclidean Yang-Mills-Higgs theory destabilises the instanton solution against 
collapse: this theory does not admit finite action solutions with non-trivial scalar. The absence 
of exact solutions does not, however, imply that one cannot make use of approximate solutions 
of the field equations. These solutions are obtained by expanding in small expectation values 
of the Higgs field at infinity. The equations of motion for the four dimensional Euclidean 
theory are 

DiF** = [$,&$], (11) 

L> 2 $ = 0. (12) 

To lowest order in the Higgs expectation value the right hand side of the equation (11) can 
be set to zero and the solution is then given by (7) with A = 0. This approximate solution 
trivially lifts to a solution of 4+1 dimensional Yang-Mills-Higgs. We will refer to this lifted 
solution, perhaps slightly inadequate, as the "constrained instantonic soliton" 6 . 

Our starting point is thus the moduli space of constrained instantonic solitons. Different 
points in this set are field configurations of different potential energy which solve equation 
of motion only up to quadratic order in the Higgs expectation value. Hence, unlike the 
situation for 1/2-BPS dyons, an effective action that describes the time evolution from one 
point in moduli space to the other necessarily involves deformations which lead to a change 
in the potential energy i.e. deformations which are not zero-modes. We will refer to generic 
deformations that satisfy Gauss' law as modes, and to the subset of these which keep the 
potential fixed as zero modes. 

Time evolution in the full field theory can be approximated by geodesic motion on the set 
of constrained instantons, provided the time evolution in the full field theory is such that a 
system which is initially tangent to this set stays near the set. In order for this to be true, 
three requirements have to be met: the valley should be completely known, the bottom of the 
valley has to be almost flat and the edges surrounding the valley should be steep. The first 
condition is necessary because there is in general no mechanism that restricts the motion to 
a subset of the valley (unless, of course, there are conservation laws at work). The second 
condition, together with the third, ensure that the static force is small, so that the system 
never develops enough kinetic energy to start climbing the mountains that surround the valley. 
Note that the whole construction is completely analogous to the one for 1/2-BPS dyons, with 
the exception that the valley is in that case precisely flat. In the general case, a slightly tilted 
valley is allowed. The fact that the set of constrained instantons is a complete set follows from 



5 It easy to check that zero-modes of the instanton are zero-modes of the corresponding instantonic soliton. 

6 In the four-dimensional theory, these approximate solutions were used by Affleck [22] to systematically 
compute instanton corrections in the Higgs phase, despite the fact that an exact solution is not available. His 
method of fixing the instanton to a given size by hand is known as the "constrained instanton method" , and 
the lowest order term is obtained in the way described above. 



the fact the scalar equation (12) admits only a single solution given a self-dual (i.e. potential 
minimising) gauge-field background. 

So let us now analyse the modes of the constrained instanton. A subset of the zero-modes 
of the instantonic soliton will also be zero-modes of the constrained instantonic soliton, while 
other zero-modes will be lifted by the presence of the non-zero Higgs. However, all of these 
modes have to be taken into account since the deformation by both types leaves us on the set 
of constrained instantons. Instanton zero-modes are very well studied in the literature (for 
a review see for instance Belitsky et al. [23]). For the charge one SU(2) instanton, there are 
eight of them: four corresponding to a change in position, one associated to the change of the 
instanton size and three more having to do with the global orientation of the instanton in the 
gauge group. The associated moduli are y l , p and the three angles 8 a (a = 1, 2, 3) which we 
have not considered so far. Our time-dependent ansatz is therefore taken to be of the form 

A o=Q)T a + y l vt 3 d,(\nH)T a , 

A = ritj d j (In H)T a , (13) 

$ — T 3 

H ' 

where the parameters p, y % and 9 now depend on time. The second term in Aq can be 
deduced by using a Lorentz boost of the static solution, or alternatively one can see that it 
corresponds to the instanton translation zero-mode. The ^-dependent terms are the global 
gauge embedding zero-modes, gauge transformed so that they sit in A . When pv and p are 
both small such that their square can be neglected, one now indeed verifies that (13) satisfies 
Gauss' law: Did^Ai = by itself and DiD l Ao = by virtue of the static equation of motion. 
In general, we need (p) 2 <C 1, (pO) 2 <C 1 and (pv) 2 <C 1 for the approximation to be valid 
(more about this issue in section 3.2). 

As the translational zero-modes are relatively uninteresting for a single dyonic instanton, 
let us set y 1 = and study the dynamics of the size and group embedding parameters. We 
insert the ansatz (13) into the action (1) and integrate over the space-like slices using the basic 
integral 

'° = / d4 *PT^4 = K> ( 14 ) 



in terms of which we have 



h = [d 4 x(Di(l/H)) 2 = I2p 2 l and I 2 = fd 4 x (d p di lni/V = 16/ . (15) 

With (Foi) 2 = 3 12 p 2 + I\ 2 and (D^) 2 = v 2 I±, the resulting moduli space action is 

Moduli = 2tt 2 fdt Up 2 + p 2 (6 a ) 2 - v 2 p 2 ) , (16) 



so that 6 a G [0, 47r). The A component always leads to a kinetic energy contribution, and 
together with the time variation of Ai they build up the first two terms in the action. These 
coincide with those of the pure instanton. The third term, which is the most interesting one, 
arises because we are in the Higgs phase, and originates from the potential energy of the 



scalar 7 . Note that, in contrast to earlier work where potentials arose, we have derived it here 
simply by inserting the explicit Higgs field from (13) in the field theory action. 

As our 6 zero modes are only valid to lowest order in a , the action above is only correct 
when the angles are small. For arbitrary angles, the factor (0 a ) 2 has to be generalised to the 
SU(2) group volume. From now on, we will for simplicity restrict our attention to a single 
angle corresponding to the U(l) subgroup selected by the Higgs value at infinity. This is a 
consistent truncation because one can check that it restricts to a geodesic submanifold. For 
this particular angle the kinetic term equals p 2 2 even for finite 0. 

The equations of motion of the model (16) then consist of a conservation law for angular 
momentum, 

0p 2 = L, (17) 

where L is an arbitrary constant, together with the evolution equation for the radius p, 

-8 p + 2 p(0 2 - v 2 ) =0. (18) 

We will solve these in the general case in the next section, but there are two special solutions 
which deserve separate attention. In the first case we look at geodesies with non-zero angular 
momentum. By choosing the value of appropriately one can cancel the potential term 
completely. In this case, (17) fixes the radius to be constant under time evolution, 

p = 0, = ±v. (19) 

This is the force balance discussed in the introduction. The geodesic rotation in the direction 
chosen by the Higgs corresponds to an exact dyonic instanton solution (with the signs corre- 
sponding to the freedom in choosing s' in (7)). The electric charge of dyonic instanton arises 
just as with the Julia-Zee dyon: it is associated to kinetic energy in a compact direction (and 
corresponds to the conserved angular momentum). 

On the other hand, when the angular momentum is zero, the radius of the instanton is 
driven to zero by the presence of the potential, and to the extent to which this approximation 
is valid, one regains the shrinking behaviour of the instantonic soliton in the Higgs phase 8 . 

2.3 Dynamics of a single dyonic instanton 

The equations of motion (18) and (17) can be solved exactly, and as we will see, the single 
dyonic instanton already shows interesting dynamical behaviour, much richer than that en- 



7 The kinetic energy f(Do&) 2 of the scalar field deserves special attention, as it diverges logarithmically. 
This is very similar to the situation for sigma model lumps (see Ward [24] , Ruback [25] and Leese [26] ) . Given 
that it comes with a coefficient v 2 p 2 , one might be tempted to interpret the divergence as a constraint on 
the dynamics of p. However, this would be incorrect for several reasons. First of all, the term comes in at 
higher order in the approximation, and could receive additional contributions. More importantly, we know 
that for vanishing electric charge, p — cannot be a solution to the equations of motion of the full theory, 
as it is forbidden by Derrick's theorem. This seems to suggest that in the approximation we make, this term 
can be ignored. Although a full resolution of this puzzle is definitely still lacking, even for sigma model lumps, 
we have verified that the cut-off method of Piette and Zakrzewski [27] indeed confirms the correctness of our 
procedure. 

8 Even in the pure Yang-Mills theory the instanton will start to shrink if it is given an initial p ^ 0, since the 
equation of motion (18) reduces to p = 0. A numerical analysis of the shrinking behaviour in pure Yang-Mills 
was given by Linhart [28, 29], who showed that the moduli space approximation result p = po + pot gets 
modifications in the full theory (the parameter p becomes a function of r and the time evolution is modified 
as well). A similar analysis could presumably be done in the Higgs phase, where corrections to the 4p = —v 2 p 
evolution of the moduli space approximation are expected to be found. 





Figure 1: Typical oscillating orbit of a single dyonic instanton. The radius is shown in 
the first plot, while the second one exhibits the deviation of the potential from its initial 
value. Observe that the potential is not constant even though the solution remains com- 
pletely within the regime of validity of the moduli space approximation. The initial data 
are pit = 0) = 1 and p{t = 0) = —0.05, while the configuration is initially rotating in the 
T 3 direction with critical velocity, 9 = v. Note that it is impossible to have a critically ro- 
tating system when p 2 goes through its equilibrium value; this follows directly by inserting 
L = vp 2 and t = —to in (23). 

countered for 1/2-BPS solitons. Upon inserting the conservation law L = 9p 2 into the equation 
of motion for p we get 



4p 



P 3 



+ v 2 p = 



c 2 . 



Multiplying by p, this can be integrated once to 

L 2 

4 (p) 2 + — + v 2 p 
p 2 

Upon further substitution of p 2 = x this becomes a standard integral, 

da; 1 . / — 2v 2 x + c 2 



(20) 



(21) 



\J — v 2 x 2 + c 2 x — L 2 



V4L 



2 v 2 - c 4 



arcsm 



vVar 



2 v 2 - c 4 



(22) 



The solution is thus a harmonic oscillator for p 2 . Eliminating the constant c in favour of the 
amplitude by using (21), the expression for p is 



P 



A sin (v(t + t )) + J-^ + A 2 = \J2A sin 2 {\v{t + t )) + J 



L 2 



+ A 2 -A. 



(23) 



In the limit L 2 = we thus recover the harmonic oscillator with frequency |w, as expected 
from (20). However, we should stress that as long as 9 ^ (or in other words, as long 
as L 2 7^ 0), the oscillating solution described by (23) does not go through zero. It is thus 
fundamentally different from the pure constrained instantonic soliton: whereas the solution 
with 6 = goes through the singular point p = in moduli space, and is therefore not 
necessarily a valid orbit for arbitrary times, the oscillations found above are genuine orbits 
that remain within the validity of the moduli space approximation. A typical orbit is depicted 
in figure 1 and compared with the 1/4-BPS dyonic instanton in figure 2. 
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Figure 2: Moduli space representation of the intrinsic dynamics of a dyonic instanton. The 
x-y plane is the moduli space of the underlying instantonic soliton. The dyon is given by the 
circular orbit at fixed height of the potential, with the angular velocity being related to the 
electric charge (light curve). More kinetic energy can be added to make the system oscillate 
around the 1/4-BPS configuration (dark curve, corresponding to the data in figure 1). 

It is also worth stressing that the generic solution does not stay at a fixed value of the 
potential. This is again different from the pure dyonic instanton with 9 = v , which does 
conserve the potential and, because the potential gradient is exactly balanced by the centrifugal 
force, has all other parameters independent of time. The only quantity which is conserved in 
the general case is the angular momentum L, related to the electric charge of the system. 

3 Dynamics of two dyonic instantons 

3.1 Metric and potential on moduli space 

With the understanding of the previous section at hand, we can now start the main part of 
our program and construct the moduli space approximation for the double dyonic instanton 
system. The modes relevant for the computation of the metric are again the zero-modes of the 
two-instanton solution. The resulting metric can be obtained from the work of Osborn [14] (see 
also the work of Bruzzo et al. [30] and Bellisai et al. [31] who used the hyperKahler quotient 
construction). In addition, there will be a potential term coming from the scalar sector. We 
will treat these two parts of the calculation in turn. 

Osborn's work [14] is focused on the calculation of the two-instanton metric determinant, 
but his paper contains enough intermediate results to extract the metric itself too. Details 
can be found in appendix A.l. Let us first consider the case when the two solitons are far 
apart. In this case the interpretation of the quaternions used in the ADHM construction (see 
appendix A.l) is straightforward: the magnitudes \v±\ and ji^l are the sizes of the instantons, 
the unit normalised quaternions fi/|f 1 1 and i^/l^l correspond to the SU(2) angles, the centre 
of mass is p and the positions are given by p + r and p — t respectively. First, consider the 
limit \t\ — ► oo. The metric then only receives contributions from the first four terms of (55). 
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In vector notation, this reads 

ds 2 = 4(d|vi| d|«i| + ki| 2 (fii) 2 + d|u 2 | d|u 2 | + H 2 (^ 2 ) 2 + (dp) 2 + (dr) 2 ) . (24) 

This correctly reduces to two copies of the flat metric of the single dyonic instanton displayed 
in the first two terms of (16). 

In case the group embedding angles are aligned the potential can be obtained from the 
't Hooft ansatz for $, 

H [r — r\Y (r — r 2 ) A 

The resulting potential has been computed by Eyras et al. [16], 

V = -(pl + pl)v 2 . (26) 

This is, as we will see, also the correct potential in the large separation limit, and the equations 
decouple completely. However, the full potential for the two-instanton case is more complicated 
and depends on the relative gauge orientation. 

When taking into account the l/(distance) 2 corrections to the metric, there are contribu- 
tions from two terms in the appendix. Interestingly, there are no corrections that depend on 
the impact parameter (i.e. that contain the angles of the quaternion r) at this order. Adding 
the second term in (57) and the first two terms (squared) in (63), and also setting the centre 
of mass coordinate to a constant, we get 



ds 2 = 4((dt; 1 ) 2 + (dt; 2 ) 2 + (d|r|) 2 ) 



+ — pr( (dv 2 df 2 )|wi| 2 + (dvidvi)|t;2| 2 + 2 (vidvi)(v 2 dv 2 ) - (vxdv^ 1 - (v 2 dviY -2(viv 2 )(dvidv 2 ) 



\r\ 2 



1 ((vidv 2 )-(v 2 dvi)) 2 . (27) 



M 2 



The term in (59) which contains an epsilon term will only contribute when the vectors V\ and 
v 2 and their variations span the four-dimensional space, i.e. they contribute only to terms in 
the metric which involve all three angles. A consistent truncation is thus to set two of these 
angles identically zero, just as we did in the single instanton case in section 2.3. In this case, 
it is convenient to use the parameterisation 

vi = p! (cos e iex + sinewy) => dui =vi(0i) — +vi(0i + f)d0i, (28) 

Pi 

and similarly for v 2 . The metric can then be written as 
ds 2 = 4 ((d Pl ) 2 + (dp 2 ) 2 + \{p\ + p 2 .) ((d0) 2 + (d0) 2 ) + |(p 2 - p 2 )d0d0 

+ 4(d|r|) 2 + — (d(i?sin0)) 2 --^(cos0(p 1 dp 2 -p 2 dp 1 ) + J Rsin0d6)\ (29) 

where we have introduced the variables 

e:=i(0! + 2 ), 

0:= i(0i-0 2 ), (30) 

i?:=pip 2 . 
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For the potential, we essentially follow the work of Dorey et al. [18]. This requires writing 
(Dj$) 2 in terms of the ADHM data and extracting the relevant parts in the large distance 
approximation. The result is 

V = v\(fZ + fZ)-±I?mn*<l>). (31) 

There are of course many terms at order 1/|t| 4 which we have ignored here. 

The metric (29) admits a Killing vector d/dQ. As expected from the requirement that our 
model can be supersymmetrised, one now indeed finds that the norm of this Killing vector 
equals (up to an overall constant) the potential, 



d 2 



V 

(32) 



,2 



Because the potential gradient contracted with the Killing vector vanishes, i.e. 

there is a conserved quantity associated to this Killing vector; it is 
L = g^BX 1 * 

= \{pl + p\ - -tt^pIpI sin 2 0j + {p\ - pj)<j) - |32 cos sin 0(pip 2 - P2P1) ^ 

@^/2 2W 1 - • -/ • 

= — T + (Pi ~ PlW - i-io COS Y S111 VKPIP2 - P2P1) ■ 

v \t\ 

The effective action and the corresponding equations of motion can now be written down 
easily using the metric (29) and the potential (31). We will refrain from spelling these out, as 
they are rather lengthy, but some parts of them will be displayed later when we discuss the 
actual dynamics. 

3.2 Validity of the approximation 

Before we analyse the orbits of the two-dyon system, let us make a few comments on the 
validity of the approximation. We have already noted that we are working in the limit in 
which v and the velocities are small, but since v and are not dimensionless they have to 
be compared to typical scales in the problem. The limit in which our approximation holds is 
therefore 

( P (tH 2 «i, [pit)) 2 «i, 

(p(t)0(t)) 2 «l, (p(t)) 2 «1, (35) 

(pW0(t)) 2 «l, (p 2 /M 2 ) 2 «l. 

In these expressions, p denotes any typical size parameter in the problem. The last condition 
could in principle be relaxed if one were brave enough to include all the higher order 1/|t| 
corrections into the metric and the potential calculation. The other conditions are intrinsic to 
the moduli-space approximation and cannot easily be avoided. 
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3.3 Analytical and numerical analysis of the orbits 

We are now ready to analyse the solutions to the equations of motion for the two-instanton 
system. From sections 2.2 and 2.3 we have learnt that there are two different types of geodesies 
for the single dyonic instanton: the BPS dyonic instanton determined by 

Pi = constant , 9i = v , (36) 

and the oscillating dyonic instanton, given by the more general solution 

p 2 = A sin (y(t + t )) + B , &\ = — 

"' (37) 

In this section we will analyse the interaction and scattering processes which as asymptotic 
states have: (i) two BPS dyonic instantons (36) or (ii) two oscillating dyonic instantons (37). 
As should be clear from the formulae for the metric and potential, we are unable to give a 
completely analytic treatment of the dynamics. The main qualitative features can however be 
extracted with some help of numerical integration 9 . 

(i) two BPS dyonic instantons: 

The first check to perform is to verify that the double dyonic instanton system at rest 
remains static at any finite separation. This is required, since (36) describes a BPS object, 
and hence it satisfies the no-force condition. One indeed easily verifies that two copies of 
solution (36) for arbitrary values of the instanton radii and arbitrary initial values of an- 
gle orientation 9i (i = 1, 2) solve the full equations of motion for the metric (29) with the 
potential (31). 

When we kick the dyonic instantons towards each other, the dynamics depends crucially 
on the value of the relative orientation <fi in the gauge group. For <fi = or <fi = n, the sizes 
remain constant and the group rotations remain undisturbed (this behaviour is in contrast to 
that of ordinary 1/2-BPS dyons, which when moving feel a force due to the difference in the 
velocity dependence of the scalar and vector forces). For the intermediate value <fi = tc/2 the 
two objects behave identically, with the radii slowly increasing as they approach each other. 
In between these values, the behaviour of the radii is more complicated. 

(ii) two oscillating instantons: 

Let us first consider two oscillating instantons which (when the instantons are far apart) 
have identical p oscillations. The equation of motion for the separation r following from 
(29,31) is given by 

-8 f - — - (pip 2 ) 2 sin 2 (f) 

\ r \ 

2 / , \ 2 



M 3 



sin0(pip 2 + P2P1) +PiP2 cos0J (3c 



2 / \ 2 

+ — jg ^ cos 0(pip 2 - P2P1) + P1P2 6 sin 0J = 



9 



The computer programs used to produce these numerical solutions can be obtained from http://www. 
damtp . cam . ac . uk/user/kp229/di . html. 
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Figure 3: Solution for which 9\ ^ 62, but the radii still satisfy p\ = p2- This is only 
possible for <fi = n/2 for all times, which implies that the group rotation of the two solitons 
is synchronised. The solitons feel an attractive force. The other initial values are p\ = p2 = 1, 
T = 10, p\ = p 2 = -0.05. A = 1.118 

From this expression one notes that evolution with <f)(t) = and identical radii, pi(t) = P2{t), 
leads to a vanishing force and hence r = constant (inspection of the equation of motion for 
shows that 0(£) = can indeed be imposed for arbitrary times). So we see that even though 
two identical and synchronised oscillating dyonic instantons are not strict BPS states in the 
full theory, they behave as BPS objects in our approximation. The oscillating motion found 
in section 2.3 solves the equations of motion of the effective action even when the solitons are 
at finite separation. 

More interesting dynamics occurs when one takes the relative angle for the gauge embed- 
ding of the instantons to be nonzero ($i 7^ 9 2 , i.e. 0^0) while keeping the p amplitudes the 
same. In this case one finds an attractive force between the instantons for all initial values 
of 0. A typical orbit is depicted in figure 3. This behaviour is difficult to extract from the 
equation of motion (38) directly (due to the non-trivial oscillations of the radii). However, we 
will see in the next section how an effective model can be derived that predicts this behaviour. 

Another, more complicated set of geodesies appears in case one has two asymptotic dyonic 
instantons which oscillate with the same p amplitudes but with opposite phases. Depending 
on the initial value of the angle one finds attractive or, in contrast to the situation before, 
repulsive behaviour. On top of this long time scale motion, there are high-frequency oscillations 
in t which were not present in the previous (equal phase) situation. The borderline case, where 
only high-frequency oscillations remain, is hard to determine, but we will see that again an 
effective model can come to the rescue. 

Finally, when the condition of equal radii is dropped, the full structure of the force between 
the solitons becomes visible. Attractive, repulsive and oscillatory behaviour of r are again 
present, but the equations are generically too difficult to analyse analytically. Typical orbits 
of this type, obtained using numerical integration, are depicted in figure 4 and 5. 

3.4 Effective action for two oscillating dyonic instantons 

Although the moduli space approach significantly simplifies the analysis of soliton dynamics 
by getting rid of high-energy modes, one can often still not do without numerical help when 
solving for the low-energy behaviour. In the case of 1/2-BPS dyons there is not much else 
one can do, but for the 1/4-BPS solitons under consideration here, we have seen that even 
in the low-energy regime there is sometimes a separation of different scales. For instance, 
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Figure 4: Generic double soliton solution of the "repulsive" type. The initial data are 
t = 6, p\ = pi = 1, p\ = —p2 = —0.05 and = 0.7. There are oscillations on top of the 
curve of r(t), but they are too small to be clearly visible in this plot. 
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Figure 5: Generic double soliton solution of the "attractive" type. The initial data are 
similar to those of 5 with the difference that now = 0.3. Again, there are oscillations on 
top of the r curve which are very small. 

the attractive and repulsive orbits discussed in the previous section exhibit high-frequency 
oscillation set by the scale of the instantons, and low-frequency dynamics set by the scale of 
separation. In the present section, we will explain how one can simplify the analysis in certain 
cases, by integrating out all but the lowest of the remaining modes. 

Let us first consider the case of two synchronised oscillating instantons, pi(t) = P2{t) with 
arbitrary and 9. As we have seen in the numerical analysis, this system always exhibits 
attractive behaviour. Strictly speaking the amplitudes of p\ and p 2 are not equal in time, 
and not even constant (the only case in which this is true is = 7r/2, given in figure 3). 
However, since the rate of amplitude decrease is much smaller than the rate of the decrease 
of the instanton separation r, one may approximate the instanton oscillation during motion 
with the free one (37). Numerical analysis also shows that the relative phase oscillates with 
increasing amplitude. However as a first approximation we will take to be a constant during 
motion. 

So we proceed by inserting the pi(t) = p^if) given by (37) and 6 = const., = const, into 
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the effective action derived from (29) and (31), and average over one period to arrive at 

f ( 2 2v 2 A 2 sin 2 <A . . 

S eS = JdtUf 2 + r J . (39) 

As expected we see that potential is attractive for all values of the angle 0. The solution to 
the equation of motion is easily found, 



^M^ . (40) 

2r 



From this expression we can see that the Higgs expectation value v does not appear in the 
separation of the scales: the time scales for p and r motion are u p = v and u T = vA/t 2 
respectively. The agreement of the low-frequency mode approximation with the full dynamics 
is excellent: when the analytical result is compared with the numerical one of figure 3, the 
errors are on the order of the width of the curve. 

It is important to stress that even though we have integrated out high-frequency behaviour 
associated to the nonzero-modes, this is by no means equivalent to ignoring them. While the 
latter corresponds to keeping the radii fixed and leads to a vanishing force (as discussed in 
part (i) of the previous section) averaging over the nonzero-modes leads to an attractive force. 

Next we consider two oscillating instantons with opposite phases 



p\ = Asin(vt) + B, t 
p\ = -A sin(i)i) + B , I 
By integrating the equation for = \{9\ — #2) one obtains 






Av 
tan(0 - (f) ) = — - cos(vt) . (42) 

Inserting (41) and (42) into the action constructed from the metric (29) and the potential (31), 
one then finds 

S= I Mr) 2 - l ^ 1 



T 2 B 2 -A 2 sm 2 (vt)l + t&n 2 (f) 

(43) 

A 2 sm 2 (vt)(l + — cos(vt)tan(<p)j - B 2 ( tan(0) - — cos(vt) J )+... 

where dots corresponds to terms which after averaging produce constant terms. By averaging 
the terms proportional to 1/r 2 one obtains again the effective force between two oscillating 
instantons. These integrations are rather tedious to do analytically, so we have restricted 
ourselves to a numerical integration for various values of the 0o parameter. One finds that the 
effective action (43) is indeed able to predict the attractive, repulsive and oscillatory phases as 
obtained in the full numerical analysis of section 3.3. Again, we should stress that integrating 
out the nonzero-modes is crucially different from simply ignoring them, and leads to much 
richer dynamics. 
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4 Discussion, conclusions and outlook 

In this paper we have exhibited several new features that arise when the moduli space method 
is extended from 1/2-BPS to 1/4-BPS soliton dynamics. The main ingredient responsible for 
these new features is the appearance of a potential, which we have explained from the point 
of view of the underlying approximately static uncharged soliton. This approach has naturally 
led us to consider dynamics for which the potential energy is not constant during evolution. 
The nonzero-mode excitations can be taken into account in a way that is consistent with 
the moduli space approximation, and provide genuinely new types of behaviour. A second 
new feature is the appearance of different scales even in the low-energy regime. This is again 
intimately related to the presence of the potential: orbits which involve motion up and down 
the potential "superposed" with motion in the flat directions exhibit such a scale separation. 
We have shown how one can integrate out the high-frequency modes among the excitations on 
moduli space, and arrive at a new effective action for the remaining modes. We were able to 
illustrate these features on the very explicit example of two dyonic instantons, as we managed 
to write down an explicit expression for the metric and the potential of this model. 

For future research, several interesting questions remain open. Firstly, the analysis of the 
dynamics is far from complete (we have only considered head-on motion of the two solitons) 
and one may want study more complicated orbits. Given the complexity of the metric, we 
have also not been able to analyse full head-on scattering that brings the solitons very close 
together. A particularly interesting question is to understand the effect of the potential on 
the 90° scattering process familiar from dynamics of other solitons. 

On a more fundamental level, it would be interesting to study the dynamical properties 
of the full supersymmetric extension of our model, as well as the quantisation of it. Previous 
studies in this direction (see for instance Bak et al. [32]) have focused on the kinematical 
properties. Finally, there are other systems for which the dynamics under influence of the 
potential can be analysed explicitly. While most of the SU(3) monopole systems have a very 
complicated metric (see Houghton and Lee [33] for details), the (1,1) system is amenable to 
direct analytical analysis. Work in this direction is in progress and will be published shortly. 
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A Appendix: technical details 

A.l Two-instanton metric 

The crucial ingredients for the computation of the metric on the two-instanton moduli space 
can be found in Osborn [14]. We will here summarise the relevant details in a coordinate 
system that is more suitable for our purposes; the reader is referred to [14] for a detailed 
explanation of the formalism. 

The square root A of the ADHM projection matrix will be parametrised as 

A = a + bx (44) 

with 



(45) 



The normalisation parameter v appearing in tfb = u± is therefore equal to one. The projector 
Poo becomes 

/l 0\ 

P^ = lim (1 - bu- 1 ^) =000 so that (1 + P^) = [ 1 ] . (46) 

\0 0/ 

We are now ready to solve the non-linear ADHM constraint a) a = fit. It leads to 

ViV 2 - v 2 v l = 2(ctt - fa) , (47) 

which is necessary to make the off-diagonal components of \x vanish (they are manifestly not 
proportional to the unit matrix in SU(2) space). This allows us to solve for a in terms of the 
other moduli: 

v 2 v x -viv 2 A 

a = T — JTT2 — + aT ='■ T (jri5 + a • 48 

4|rp 4|t| z 

We will set a = in order to be able to use the interpretation of the moduli as given by Dorey 
et al. [18]. In order to make contact with the work of Osborn [14], we have to find a real 
function A such that a above satisfies 

v 2 vi -2fa = 2A|r| 2 . (49) 

One finds that 

V 2 V X + v x v 2 
X= 2|r|* • (50) 

We will see in a moment that this is proportional to 5 a b in spinor space. It is also useful to 
have a compact expression for A, 

A = 4«a mn = -A. (51) 
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The metric for the zero-modes can be obtained from (4.14) of Osborn [14], 

(C", C) = (S'a, 5a) - Ak'kN^ 1 , (52) 

:={C',C)i :=(C',C) 2 

where the inner product of the first term is given by 

(S'a, 5a) = tr (S'a^(l + P^Sa) . (53) 

The first term in (52) leads to a diagonal contribution to the metric, which is easy to see by 
computing 

5'a ] (l + P OQ )5a = 

'2 5'v 1 5v 1 + (S'p + S'f) (Sp + St) + 5'a5a ■■■ \ 

. (54) 
2S'v 2 Sv 2 + (S'p-S'f)(Sp-Sr) + S'a8crJ 

The trace then eliminates the SpSr cross-term, and the result for the first term in the metric 

is 

(C, C) 1 = 2 tr 2 (S'v 1 Sv 1 + 5'v 2 5v 2 + S'pSp + SfSr + 5'aSa) . (55) 

The final step consists of eliminating the a in the last term. The variations of a and a are 

S'a = S'Af-^—r + AS'f -^ - KfS'irf ' 



4|r| 2 4|r| 2 v 7 4|r| 4 ' 

(56) 
11 1 

5a = 5tA—-t + tSA—^ - tA5(tt)— — . 
4|t| 2 4|r| 2 v 4|r| 4 

The product thus leads to nine different terms for which we have to work out the conversion 
of quaternions to vectors. 

(S'a) (5a) = S'AfSfA-^—- + 5' A5 A— ]— - S'AAS(rf) ' 



16|r| 4 16|r| 2 v y 16|t| 4 

+ AS'f St A— — + AS'f tSA— — - AS'frAS(rf) —— & (57) 

lo|r| 4 lo|r| 4 lo|r|° 

- AfS'(rf )StA-^- e - AS'(rf)SA-^- + AAS'(rf )S(rf )-J- |? . 
1d|t| d lo|r| 4 loir) 

In order to limit the number of a matrices appearing in the trace, it is convenient to rewrite 
S(rf) right from the beginning in vector notation as 

S(rf) = 21 2 (r • St) . (58) 

In the main text we ignore all 1/|t| 4 terms, so there is only one term remaining, 

| ti(S'ASA) = 

(S'v 2 ■ Sv 2 )\vi\ 2 + (vi ■ S'vi)(v 2 ■ Sv 2 ) + (vi ■ Svi)(v 2 ■ S'v 2 ) + (S'vi ■ 5vi)\v 2 \ 2 

- (S'v 2 ■ vi)(Sv 2 ■ vi) - (S'vi ■ v 2 )(Sv 1 ■ v 2 ) - (v 2 ■ vi)(S'vi ■ Sv 2 ) - (v 2 ■ v 1 )(Sv 1 ■ S'v 2 ) 

- (v™8'vl8v k 2 v[ + v k 2 Sv[S'v^)e mnkl . 

(59) 
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Now we focus on the second part, (C, C) 2 . For this we have to compute 



a ] 5a - {a [ 5a) T = K = t 2 ® ( 



k 
-k 



(60) 



For the matrix element K 12 one finds 

K12 = ViSv 2 — v 2 8v\ + 2f5a — 2a5r . (61) 

Inserting the solution (48) of the ADHM constraint and using (58), this becomes 

K 12 = \{v\5v 2 + 8v 2 v 1 - v 2 5vi - Sviv 2 ) - |(r-f^^ r + ^ fT r~^) ■ (62) 

As expected this is indeed proportional to the identity matrix in 577(2): the terms that only 
have Vi are proportional to o m <y n + a n a m , while the other ones are proportional to d'mnO'ki, 
both of which (by (69d) and (69a) respectively) are proportional to the identity. Rewriting k 
in vector notation, we find 

k = (ui • 5v 2 ) - (v 2 ■ Svi) - 7^(e mraM «r fc (5r' - (v 2 ■ r)(«i • St) - {v x ■ t)(v 2 ■ Srfj . (63) 

From this we now have to compute k'k(NA)~ 1 - This requires the matrix Na, which is given 
in (3.30) of Osborn's paper, 

N A = M 2 + H 2 + 4(|r| 2 + \a\ 2 ) = \v x \ 2 + \v 2 \ 2 + 4|r| 2 + \v 2 v x - v 2 v x \ 2 . (64) 

The result for k'kN^ 1 can now easily be read off. 

A. 2 Conventions and sigma algebra 

This appendix collects some useful expressions for Euclidean sigma matrices (indices m,n, . . . 
take values 0, 1, 2, 3 but the zeroth direction has positive signature too). The basic definitions 
in terms of the Pauli matrices r l are 

a m := (1, ir*) , (65a) 

a m := (1, -ir l ) , (65b) 



a 



a 



mn If ^.rn — n ^n^m 

■~ 4 

mn . 



(a m a n - a n a m ) , (65c) 



i(aV-aV), (65d) 



1 



a 



mnp — l(m-np _ p-n m 



2 



(a m a n a p - a p a n a m ) , (65e) 



a 



mnp ._ 1 l-mn-p _ -p n-m 



2 



(a m a n a p - a p a n a m ) . (65f) 



The sigma matrices satisfy the algebra 

aia m = 2a Lm + 5 lm , (66a) 

o\(j m = 2 aim + Si m ■ (66b) 
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O'mn 


— 


1 
2^-mnpqOpq i 


0~mn 


— 


2~£mnpq&pq j 


•mnp 


— 


rnnpq 

e o q , 


— mnp 


— 


rnnpq zf 

c Oq . 



and one can derive analogues for products of more sigma matrices. 

G- n a kl = \a nkl + <5 n V] = \e nkl ic q + S n[k a l} , (67a) 

"mynkl ^mnkl ^nklq"mq j ^0 1 DJ 

CmnO'kl = —^SmkSnl + l&ml&nk + ^e" 1 ™ , (67c) 

0~mn<3~kl = —^mk^nl + l&ml&nk ~ -e™ - ™ . (67d) 

The epsilon symbols (which satisfy £0123 = 1) arise because of the following duality flip rela- 
tions, 

(68a) 
(68b) 
(68c) 
(68d) 
Finally, one can of course trace these expressions to arrive at 

tr 2 (a mnp ) = 2e mnp0 , (69a) 

t r 2 [O'mnO'kl) = ~ ^nik^nl + ^ml^nk + 2~ e "" ' (69b) 

tf2 {O'mnO'kl) = _ ^mk^nl + ^ml^nk ~ 2^ ' (69c) 

tr 2 {a m a n + a n a m ) = A5 mn . (69d) 
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